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Abstract 

First measurements of cross sections for isolated prompt photon production in 
deep inelastic ep scattering have been made using the ZEUS detector at the 
HERA electron-proton collider using an integrated luminosity of 121 pb~^. A 
signal for isolated photons in the transverse energy and rapidity ranges 5 < 
< 10 GeV and —0.7 < rj'^ < 0.9 was observed for virtualities of the exchanged 
photon of > 35GeV^. Cross sections are presented for inclusive prompt 
photons and for those accompanied by a single jet in the range E}^^ > 6 GeV and 
— 1.5 < 77-''^* < 1.8. Calculations at order a^Ug describe the data reasonably well. 
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1 Introduction 



Isolated photons in the final state with high transverse momenta are a direct probe of the 
dynamics of hard subprocesses in high energy collisions, since these 'prompt' photons are 
largely insensitive to the effects of hadronisation. Prompt photons have been studied in a 
number of hadronic experiments. Early evidence for such processes came from the R806 
experiment at the CERN ISR [1]. More recently, the CDF and D0 experiments at the 
Tevatron collider have performed a number of QCD tests using prompt photons [2-7]. In 
previous ZEUS publications, the production of prompt photons in photoproduction has 
been studied [8-10]. In the present letter, for the first time, prompt photon measurements 
in deep inelastic scattering (DIS) are reported, both inclusively and accompanied by jets. 
These processes test QCD in a new way by studying processes containing two different 
hard scales, Q^, the exchanged photon virtuality, and E^, the transverse energy of the 
emitted prompt photon. 

Prompt photons are produced in DIS at lowest order in QCD, as shown in Fig. 1. These 
processes have been calculated to order 0{a^as) by Gehrmann-DeRidder, Kramer and 
Spiesberger [11], including interference terms for initial- and final-state radiation from the 
electron. In contrast, leading- logarithm parton-shower Monte Carlo (MC) models do not 
naturally predict events with two hard scales. 

In this letter, results are presented for the process ep e'-fX, where X is anything, and 
for ep —>■ 67 -|- jet -|- Y, where Y does not contain further jets within the acceptance of the 
measurement. Comparisons are made to MC predictions and also to 0{a^as) calculations 
for the photon-jet final state. 

2 Experimental set-up and event selection 

A data sample corresponding to an integrated luminosity of 121 pb~^ was used, taken 
between 1996 and 2000. This sample is the sum of 38 pb~^ of e~^p data taken at a centre- 
of-mass energy of 300 GeV and 68 pb"^ taken at 318 GeV, plus 16 pb~^ of e~p data taken 
at 318 GeV. A single set of results is presented for this combined sample. The MC cross 
sections (see Section 3) differ by under 4% at the two centre-of-mass energies, well within 
the precision of these measurements. Differences between the cross-sections for e~^p and 
e~p collisions are expected to be negligible [12]. 

A description of the ZEUS detector is given elsewhere [13]. Of particular importance in 
the present work are the uranium calorimeter (CAL) and the central tracking detector 
(CTD). 
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The CAL [14] has an angular coverage of 99.7% of Air and is divided into three parts 
(FCAL, BCAL, RCAL), covering the angular ranges 2.6° - 36.7°, 36.7° - 129.1° and 
129.1° — 176.2°, respectively^. Each part consists of towers longitudinally subdivided into 
electromagnetic (EMC) and hadronic (HAC) cells. The electromagnetic section of the 
BCAL (BEMC) consists of cells of 23.3 cm length azimuthally, representing 1/32 of the 
full 360°, and width of 4.9 cm in the Z direction at its inner face, at a radius of 123.2 cm 
from the beam line. These cells have a projective geometry as viewed from the interaction 
point. The profile of the electromagnetic signals observed in clusters of cells in the BEMC 
discriminates between those originating from photons or electrons^ and those originating 
from neutral-meson decays. The CAL energy resolutions, as measured under test-beam 
conditions, are a{E)/E = 0.18/v^ for electromagnetic showers and a{E)/E = 0.35/v^ 
for hadrons, with E in GeV. 

The CTD [15] is a cylindrical drift chamber situated inside a superconducting solenoid. 
Using the tracking information from the CTD, the vertex of an event can be reconstructed 
with a resolution of 0.4 cm in Z and 0.1 cm in X,Y. In this analysis, the CTD tracks 
are used to reconstruct the event vertex, and are also used in the selection criteria for 
high-ET photons. 

The luminosity was determined from the rate of the bremsstrahlung process ep — ^ e'jp, 
where the high-energy photon was measured in a lead-scintillator calorimeter [16] located 
at Z = —107m. 

The DIS events were selected online using a trigger based on energy deposits in the CAL 
consistent with a scattered electron. Offline, events which passed DIS cuts similar to 
those used in previous analyses [17] were selected. In addition a photon candidate was 
required. The value of Q^, as reconstructed from the final-state electron, was required to 
be above 35 GeV^. The energy of the scattered electron was required to be above 10 GeV 
and its polar angle in the range 139.8° to 171.9°, in order to be well measured in the 
RCAL and well separated from the photon candidate. Events were required to have a 
reconstructed vertex position within the range \Z\ < 40 cm and 35 < 5 < 65 GeV, where 
6 = ~ cos6'j), Ei is the energy of the ith CAL cell, 6i is its polar angle and the 

sum runs over all cells. 

For the subset of events used in the photon-jet study, jets were reconstructed from CAL 
cells using a cone algorithm with radius 0.7 [18] in the laboratory frame. Corrections for 
energy losses, principally due to uninstrumented material in front of the CAL, were eval- 

^ The ZEUS coordinate system is a right-handed Cartesian system with the Z axis pointing in the proton 
beam direction, refered to as the 'forward direction', and the X axis pointing left towards the centre 
of HERA. The coordinate origin is at the nominal interaction point. 

^ Hereafter 'electron' refers both to electrons and positrons unless specified. 
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uated using MC simulated events, and were typically +(10-15)% for jets with measured 
energy above 6 GeV [10]. 

3 Monte Carlo event simulation 

The MC programs Pythia 6.206 [19] and Herwig 6.1 [20] were used to simulate prompt 
photon emission for the study of event-reconstruction efficiency. In both generators, the 
partonic processes are simulated using leading-order matrix elements, with the inclusion 
of initial- and final-state parton showers. Fragmentation into hadrons is performed using 
the Lund string model [21] in the case of Pythia, and a cluster model [22] in the case of 
Herwig. The events generated using the Pythia and Herwig programs were used to 
correct for detector and acceptance effects. The corrections provided by Pythia were used 
as default and those given by Herwig were used to estimate the systematic uncertainties 
due to the treatment of the event dynamics and of parton showering and hadronisation. 
The detector response to photons and neutral mesons (7r° and rj) was simulated by using 
single-particle MC generated events. 

The generated events were passed through the ZEUS detector and trigger simulation 
programs based on Geant 3.13 [23]. They were reconstructed and analysed by the same 
programs as the data. The jet search was performed using the energy measured in the 
CAL cells in the same way as for the data. The same jet algorithm was also apphed to 
the final-state particles. 

To study the effects of electron radiation, simulations were made of deep inelastic scat- 
tering events using the Heracles 4.6.1 [24] program with the Djangoh [25] interface 
to the MC generators that provide the hadronisation. The coUinear radiative corrections 
were found to be small in the kinematic region of this analysis and were neglected. 

4 Photon candidate selection 

The identification of events containing an isolated prompt photon candidate follows closely 
the approach used in previous analyses [8-10]. Events were selected on the basis of an 
isolated photon candidate detected in the BCAL. The algorithm selected predominantly 
electromagnetic clusters of cells within a small angular cone. Initially, larger electromag- 
netic clusters than are typical of a single photon were accepted to estimate backgrounds. 
Use of shower shapes as a discriminant, as described below, allowed subtraction of the 
backgrounds due to 7r° and rj production. 
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It was required that the reconstructed transverse energy of the cluster satisfied Erj. > 
5 GeV and the pseudorapidity satisfied —0.7 < rp < 0.9. The cut E'^ < 10 GeV was 
imposed to ensure that the 7r° and rj subtraction method was effective. 

The photon candidate was well separated from the scattered electron. Monte Carlo sim- 
ulations and 0{a^as) calculations (see Section 6.2) show that for electrons in the range 
defined in Section 2, most photons radiated from the electron fall outside the prompt- 
photon acceptance used in this analysis, though they still give an important contribution 
to the cross section in the kinematic region of the measurement. 

To reduce backgrounds, the photon-candidate cluster was required to be isolated by de- 
manding Ar > 0.2, where Ar = ^A^^ + Arj"^, the distance to the nearest reconstructed 
track in 77 — space. It was further required that E]^/E^^'^ > 0.9, where E^'^'^ is the energy 
within a cone in 77 — of radius 1.0 around the photon candidate. This energy isolation 
requirement suppresses the contribution from photon candidates produced within jets. 
Deeply virtual Compton scattering (DVCS) events were removed by demanding at least 
two tracks reconstructed in the CTD, since in DVCS the final state seen in the detector 
consists only of a photon and an electron which are well separated [26,27]. 

The selected candidates were still dominated by neutral mesons, such as 7r° and 77, which 
decay to photons. The single-photon signal was statistically extracted from the back- 
ground using BEMC energy-cluster shapes. The first distribution considered was that 
of {SZ), where {6Z) = S(i?ceii|^ceii — Z\)/TiEce\\- Here £^ceii is the energy deposited in 
a BEMC cell, Zceii is the cell number measured in the Z direction and Z is the energy- 
weighted mean of Zcew- Figure 2a shows the {6Z) distribution for data, together with a 
fit based on photon shower shapes and a simulation of single particles in the detector (vr" 
and 77). Clear peaks are visible at (SZ) ~ 0.15 due to single photons and (SZ) ^ 0.5 due 
to vr" — > 77, as well as a tail due to the decays of heavier particles to two or more photons. 

The photon shower shapes used were derived in two ways: from DVCS data [27], and 
from single-photon MC simulation. In Fig. 2, photons found in DVCS data events are 
shown. The results of the two shower-shape methods gave indistinguishable background 
subtractions and differed only by an overall scale factor of 5% on the acceptance of the 
prompt-photon signal. The DVCS method gave the higher acceptance, as the DVCS 
single-photon showers are slightly narrower than those from the MC showers. The MC 
method was used in this analysis, because of the higher statistics available. This allows 
rapidity and energy dependences of shower shapes to be modelled; a scale correction of 
5% was then applied. 

The ?7 contribution was determined from a fit to the {SZ) distribution above 0.65. After 
removing candidates with (SZ) > 0.65, the final background subtraction was performed 
using the variable /max, defined as the ratio of the energy of the highest-energy cell in an 
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electromagnetic cluster to the total cluster energy. When incident on the BEMC, single 
photons form narrow clusters, with most of the energy going into only one cell, giving 
an /max distribution peaked close to unity. Because of the projective geometry of the 
BEMC, a photon entering at the boundary between two cells typically has /max — 0.5. 
Thus the /max distribution for single photons peaks close to 1.0 and extends down to 0.5. 
In contrast, the neutral mesons decay to more than one photon, forming larger clusters 
in the BEMC. In each bin of a plotted physical quantity, events were divided into two 
classes, with high and low values of /max respectively. From the number of events in each 
class, as well as the ratios of the corresponding numbers for the /max distributions of the 
single-particle samples, the number of events in the given bin was evaluated [8]. 

A total of 1875 events with {bZ) < 0.65 were selected, of which 877 have /max > 0.75, 
yielding a signal of 572 and a background of 1303 events. The fits and signal extraction 
procedure were repeated for each bin of each distribution. 

Studies based on single-particle MC samples showed that the photon energy measured 
in the BCAL was on average less than the true value, owing to energy loss in the unin- 
strumented material in front of the BCAL. To compensate for this effect, a correction of 
typically 0.2 GeV was added to the photon energy [10]. 



5 Systematic uncertainties 

The following sources of systematic uncertainty were investigated: variations of the nom- 
inal /max spectra for the photon affecting the signal extraction; change in the detector 
energy scale calibration by ±3%, reflecting the overall energy scale uncertainty; and a 
change in the energy cut in both MC events and data by ±10% for photons. This last un- 
certainty is motivated by the r.m.s. differences between hadron-level generated and recon- 
structed energies. Also included as a systematic uncertainty is the difference in estimated 
acceptance between Herwig and Pythia, which is mostly well below the statistical un- 
certainty. A change of ±20% in the hadronic energy cut for photon-jet events for both 
data and reconstructed Monte Carlo events, representing the r.m.s. difference between 
hadron-level and reconstructed jet energies was considered as an additional systematic 
uncertainty. The uncertainty of 2.2% on the luminosity measurement was neglected in 
the differential cross sections but included in the total cross sections. 

The method used for background subtraction is more sensitive to the shape of the /max 
distribution of the background than to that of the signal. The background shape is 
relatively insensitive to the tt^/t] ratio and hence the results using DVCS and MC photons 
are very similar. A study was made of the effect on the results of the fact that the fits fall 
below the data at high (SZ). This is due to events with large E^, where the contribution 
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of events with more than one 7r° with a multi-7r° invariant mass above the rj mass is hkely 
to be important. A fit was made to the high-E'^ data excluding the region {6Z) > 1.0. 
The change in the extracted signal was well below the statistical uncertainty in the bin. 

6 Results 

6.1 Inclusive prompt photon production 

The cross section for inclusive prompt photon production, ep — > e'-fX, has been measured 
in the following kinematic region: > 35 GeV^ > 10 GeV, 139.8° < 6^ < 171.8°, 
—0.7 < 7]'^ < 0.9 and 5 < E'J, < 10 GeV, with photon isolation such that at least 
90% of the energy found in an 77 — cone of radius 1.0 around the photon is associated 
with the photon. The measured cross section is 

a{ep e-fX) = 5.64 ± 0.58(stat.)+[|;72(syst.) pb. 

The predicted cross sections from Pythia and Herwig are lower than the data by factors 
of approximately 2 and 8, respectively. Figures 3a and 3b show the measured rapidity 
and transverse energy distributions, compared to MC predictions normalised to the data. 
The data are also presented in Table 1. Both Pythia and Herwig describe the 
spectrum and Herwig describes the rapidity well. Figure 3c shows the distribution 
of the data, again compared to MC predictions. The agreement of Pythia with the data 
is reasonable, but Herwig fails to describe the measured spectrum. As discussed 
in Section 6.2, the 0{a^as) calculations suggest that the discrepancies between Pythia 
and the data in the rate and photon rapidity distribution may be due to the fact that 
wide-angle initial- and final-state radiation from the electron are not included in the MC 
calculations. 

6.2 Prompt photon plus jet production 

Owing to divergences in cross-section calculations for prompt photons, a comparison to 
0{a^as) QCD predictions in DIS can be made only when there is a jet accompanying an 
isolated prompt photon. Jets were reconstructed as described in Section 2. For events 
satisfying the criteria for isolated prompt photons described above, jets were counted only 
if they had E^^^ > 6 GeV and —1.5 < 77^'^* < 1.8. The measured total cross section for 
photon plus a single jet within this kinematic region is 

a{ep 67 + jet + Y) = 0.86 ± 0.14(stat.):^J:34(syst.) pb. 



6 



Figure 4 shows the differential cross sections for 'prompt photon plus one jet' events, 
together with MC predictions. The data are also presented in Table 2. The transverse 
energies of the photon and the jet are well described by the MC calculations. Herwig de- 
scribes the photon rapidity well but the jet pseudorapidity peaks at lower values. Pythia 
describes the jet pseudorapidity well, but the photon rapidity peaks too far forward, as 
was also the case for inclusive photons. 

Figure 5 shows the same data as Fig. 4, compared to the 0{a^as) parton-level calculations 
of Kramer and Spiesberger [12]. These include all possible initial- and final-state single 
photon and gluon radiation, together with appropriate vertex corrections, and their inter- 
ference terms. Higher-order effects, such as coUinear bremsstrahlung in the same event as 
a hard non-collinear photon, estimated to be a 4% effect, are omitted. These calculations 
use the phase-space-slicing method to cancel the infrared and coUinear singularities. The 
MRST parton distributions [28] were used for the parametrisation of the proton struc- 
ture. Parton-to-photon fragmentation functions were taken from Bourhis, Fontannaz and 
Guillet [29]. The renormalisation scale was chosen to be the transverse energy of the 
jet. The effect of changing this scale up or down by a factor of two, to estimate the 
possible contribution of unknown higher-order terms, is shown in Fig. 5. The predicted 
total cross section for the mixture of energies and beam charges used in this analysis is 
1.33 ± 0.07 pb, where the uncertainty corresponds to the change in the result when 
the renormalisation scale is varied by a factor of two. This parton-level calculation is 
compatible with the data. 

By definition, the 0{a^as) parton-level calculation does not include the effects of hadro- 
nisation. Hadronisation effects were investigated by comparing the parton-level Pythia 
and Herwig distributions with the hadron level. The effect of hadronisation would be 
to reduce the predictions by 30% to 40%. Because of the overall poor description of the 
data by the MC simulations, hadronisation corrections were not applied to the 0{a^as) 
calculation. 

The 0{a^as) calculation shows that 65% of photons are emitted by the electron, con- 
centrated at low photon rapidities, and the rest by quarks. The photon rapidity and jet 
pseudorapidity distributions for the latter component resemble the Pythia predictions, 
which include only such photons. Interference between these processes contributes only 
2% to the total. The transverse-energy distributions of the two processes are similar. The 
0{a^as) calculation predicts a higher jet cross section at forward pseudorapidity and at 
low EZ than is seen in the data. 
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7 Conclusions 



The first observation of prompt photon production in deep inelastic scattering has been 
presented, together with distributions for accompanying jets. Leading-logarithm parton- 
shower Monte Carlo models for photon emission by quarks (Pythia and Herwig) are 
each able to describe some but not all of the features of the data. Both describe the 
transverse energy distribution well and Herwig describes the photon rapidity well. Both 
models predict too low a cross section. 

The results have been compared to an 0{a^as) parton- level calculation for ep —>■ e7+jet + 
Y in the acceptance region of this measurement. The level of agreement is satisfactory 
in photon rapidity and jet transverse energy but only fair for photon transverse energy 
and jet pseudorapidity. The total predicted cross section is consistent with the measured 
value. 
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^7 


0.(7 j GLTj l^puj 


n 7 tn n 


A Qc; -un 78+0-51 

4. JO itU. ^8_o.63 


-0.3 to 0.1 


5.20 ±0.75l^i^ 


0.1 to 0.5 


2.12 ±0.69l[j:i 


0.5 to 0.9 


1.82 ±0.66l^:^^ 


(GeV) 


da/dE^ (pb GeV^i) 


5.0 to 6.0 


1.25 ±0.40+°:°^ 


6.0 to 8.0 


1.40 ±0.20l[J;[i^ 


8.0 to 10.0 


0.79 ±o.iol[!:J? 



Table 1: Differential cross sections for inclusive production of isolated photons with 
—0.7 < rp < 0.9, for 5 < E'^ < 10 GeV . The first uncertainty is statistical, the second 
is systematic. 





dcr/dr]'^ (pb) 




dfj/dr/j^* (pb) 


-0.7 to -0.3 
-0.3 to 0.1 
0.1 to 0.5 
0.5 to 0.9 


0.81 ±0.20+!^:^? 
0.77 ±0.19+!^:^? 

0.16 ±o.i6+!]:i? 

0.30 ±0.15toil 


-1.5 to -0.84 
-0.84 to -0.18 
-0.18 to 0.48 
0.48 to 1.14 
1.14 to 1.8 


0.087 ±0.043+[J;[J^^ 
0.118 ±0.068l|];?[J 
0.47 ±0.13+!^:?^ 

0.24 ±o.io+!^:^S? 

0.41 ±0.12t°:i° 


(GeV) 


da/dE^ (pb GeV^i) 


E^' (GeV) 


da/dE^^' (pb GeV"i) 


5.0 to 6.0 
6.0 to 8.0 
8.0 to 10.0 


0.136 ±omotoM3 

0.258 ±0.050;[J;^^? 
0.107 ±0.029;[5;g^? 


6.0 to 8.0 
8.0 to 10.0 
10.0 to 12.0 
12.0 to 16.0 


0.177 ±0.044+°:°^^ 
0.132 ±0.049l[!;!J^^ 
0.082 ±0.032t|];°?^ 
0.046 ±0.0311°:°^? 



Table 2: Differential cross sections for production of isolated photons plus one jet 
with -0.7 < r/T < 0.9, for 5 < < 10 GeV and -1.5 < rf''^ < 1.8, for > 6 GeV. 
The first uncertainty is statistical, the second is systematic. 
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Figure 1: The lowest-order tree-level diagrams for prompt photon production in ep 
scattering. Vertex corrections enter at the same order. 
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Figure 2: (a) Distribution of {6Z) for prompt photon candidates in selected events, 
(h) Distribution of f^i^i after a cut on {6Z) < 0.65. Also given are fitted distributions for 
Monte Carlo rj mesons, tt^ + rj and vr*^ + ?] + 7 (where the 7 is taken from DVCS data), 
with similar selection criteria and Ej, spectrum to the observed candidates. 
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Figure 3: Inclusive prompt-photon differential cross section (a) in rapidity, (h) in 
transverse energy, in the range —0.7 < rf < 0.9 and 5 < Ej, < 10 GeV . The inner error 
bars are statistical while the outer represent systematic uncertainties added in quadrature, 
(c) Distribution of Q'^ . In each case the histograms show MC predictions, normalised to 
data. 
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Figure 4: Cross section for prompt-photon-plus-jet production differential in (a) photon 
rapidity, (h) photon transverse energy, (c) jet pseudorapidity, (d) jet transverse energy, 
for events with a photon in the range —0.7 < rf < 0.9 and 5 < Ej, < 10 GeV and one 
jet in the range —1.5 < if'^^ < 1.8 and Ej!^ > 6 GeV . The inner error bars are statistical 
and the outer represent systematic uncertainties added in quadrature. The band around 
the data points shows the effect of calorimeter energy-scale uncertainty. The histograms 
show Monte Carlo predictions, normalised to the data. 
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Figure 5: Cross section for prompt-photon-plus-jet production differential in (a) photon 
rapidity, (b) photon transverse energy, (c) jet pseudorapidity, (d) jet transverse energy, 
for events with a photon in the range —0.7 < rf < 0.9 and 5 < E'j, < 10 GeV and 
one jet in the range —1.5 < rf'^^ < 1.8 and E'!^' > Q GeV . The inner error bars are 
statistical while the outer represent systematic uncertainties added in quadrature. The 
band around the data points shows the effect of calorimeter energy scale uncertainty. The 
boxed band shows the parton-level predictions of Kramer and Spiesberger including the 
effect of renormalisation scale uncertainty. The single line indicates their prediction of 
the contribution of photons radiated from the quark line. 
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